home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CSIFA.z / CSIFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCSSSSIIIIFFFFAAAA((((3333FFFF))))                                                            CCCCSSSSIIIIFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CSIFA   - CSIFA factors a complex symmetric matrix by elimination with
  10.      symmetric pivoting.
  11.  
  12.      To solve  A*X = B , follow CSIFA by CSISL.  To compute  INVERSE(A)*C ,
  13.      follow CSIFA by CSISL.  To compute  DETERMINANT(A) , follow CSIFA by
  14.      CSIDI.  To compute  INVERSE(A) , follow CSIFA by CSIDI.
  15.  
  16.  
  17. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  18.       SUBROUTINE CSIFA(A,LDA,N,KPVT,INFO)
  19.  
  20. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  21.      On Entry
  22.  
  23.      AAAA COMPLEX(LDA,N)
  24.         the symmetric matrix to be factored.
  25.         Only the diagonal and upper triangle are used.
  26.  
  27.      LLLLDDDDAAAA INTEGER
  28.         the leading dimension of the array  A .
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAA a block diagonal matrix and the multipliers which
  34.         were used to obtain it.
  35.         The factorization can be written  A = U*D*TRANS(U)
  36.         where  U  is a product of permutation and unit
  37.         upper triangular matrices , TRANS(U) is the
  38.         transpose of  U , and  D  is block diagonal
  39.         with 1 by 1 and 2 by 2 blocks.  KVPT    INTEGER(N)
  40.         an integer vector of pivot indices.
  41.  
  42.      IIIINNNNFFFFOOOO INTEGER
  43.         = 0  normal value.
  44.         = K  if the K-th pivot block is singular.  This is
  45.         not an error condition for this subroutine,
  46.         but it does indicate that CSISL or CSIDI may
  47.         divide by zero if called.  LINPACK.  This version dated 08/14/78 .
  48.      James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.  Subroutines and
  49.      Functions BLAS CAXPY,CSWAP,ICAMAX Fortran ABS,AIMAG,AMAX1,REAL,SQRT
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.